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SYMBOLIC INTEGRATION - II 


by Joel Moses 


In Cnis memo we describe the current state of the inte- 
gt-ation program originally described in AX Memo 9? (MAC-K-310). 
Familiarity with Memo 97 ia assumed. Some ot the algorithms 
described in that memo have been emended. Certain new algorithm- 
and a s-mple integration by.parts routine have been added. The 
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X. x & Elem(x < ), where c, V ^ are integers, and where 
k « g.c.d. of the k- divides (c + 1). 

Substitute y * x.* 

) x"/(x ,A + i)dx becomes 1/4 ;y/(y 4 1 ) dy 

M/x ~ y dx, becomes 1/2 \ Vy -T dy 


II. (Rational) R(x). Chains to a program written by Manove 

which integrates rational functions by a method described in 
Hardy ? 


111. Elem(base 
constants. 


a rb .X a , f ijX 

. , base^ , ...)> where base^, a *, b', are 


a. base: *■ ^ is converted to base/ * * ^ J W* l 




b. y ■ base, is substituted 


r 


J 


e^(2 4 3e^ )dx becomes jl/(2 + 3y A )dy 
- *#■ 
ic/e^dx becomes ' y**® 4 dy 


IV* 


V* f a 

(Chebyschev) Ax (c^ 4 c^x* ). A,r,r\j ,q,Cj ,c^ are constants. 
Convert A to A/q, let r. * r 4- 1 — 1 

1 q 

a. integer, > 0 
Substitute y = c^~ ”+~ e A $ 

b. rj integer, r. rational number with denominator d i 

1 a *U. 

Substitute y * x*^- 


c. r« integer, r, < 0, r. rational with denominator d 

In 

Substitute y * (c^ 4 c ^ k *)"** 

d. rj + r a is an integer ^ . 

Substitute y ■ /c^ 4 c^x^ '> /Jj. 

I - rr y 



v. 


F.leffi(x, ^ax + b 


\ c y. ■+ o 


ex 4 b 
CX -4' ft 
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I-et k fc least cor.non multiple .it me n ( * 

Substitute, y « -V ax -f I 

V ex 4- o 


jcof Va dx becomes 2y 

i J 

f x "/x"" + 1 dx be "ones 


cos ;v<dv 
* *■ • 

j- ^ - i)y" l -y 


vi. r (x,1 a + bx J 

a. a > 0, b > 0 Substitute car. (y ) c Va/b x 
j 1/(4 + x^Jdx becomes j^l/2 sec(y)dy 
’ b. a > 0, b < C Substitute sin(y) « Va/b x 
fx/(l - x J )^dx becomes ^fsin(y)dy 

c. a < 0, b> 0 Substitute sec (y) *'V-a/b x 

VII. R(x, 'Vex* + bx 4 a) 

a. c > 0 Substitute 

b. a ^ 0 Substitute 


y = 'Vc" x +^ cx + bx + a 
y «A/cx J «f bx + a - Va)/x 


VIII. 


R(x> 

a. 


b. 

c. 


/■ 


F(S(x)), where R(x)dx is rational 




F = log 


& /• rt-y — 

Solution is lcg(S(x)) ^R(x)dx - j | ;R(x)dx | S*(x)/S (x)cix 


r ^ f 

jX. lc-g(x)ax becomes 1/2 x' iog(x) — x/2 dx 

F c arc tan 

Solution is arctan(S(x)) 

F = arcs in r ^ , 

Solution is arcs in (S (x)) ,R(x)dx — j R(x)dx S(x)/(1 - S 4 (x)) J ax 

y _y 1/ 

jx** arcsin(x)dx becomes 1/3 x 3 arcsin(x) - J1/3 x"V(l - x A )*dx 




R(x)dx — ( R(x)dxjS(x)/(l + S (x))dx 


J U 


"J 


IX. R(x) Elem (log(a + bx)) 

Substitute y *» log (a + bx) 

jlog(x)/(l + log(x)) dx becomes jc*' y/(l 4 y) dy 



Klem(F (a + bx)) where the F: are trigonometric functions. 

t i 

I. Preparatory step y « a by yielding 1/b £lem(F-(y)) 

Note: Some code has been written for the case in vnich 
the arguments of the t rigor, om* trie functions arc linear 

but not identical. 


II. (sin (y)cos (y)dy 

J r 

a. n < n; transform into ) (1/2 sin 2>) (1/2 + i; 


.'2 c.cz 2 y) d y 


b. m2: n; transform into ^1/2 sin 2y) (1/2 - 1/2 cos 2v)'tly 

. > ^ <■ 

jeos (x)dx becomes ;(1/2 + 1/2 cos 2y)dy 

III. i) Transform all trigonometric functions into sinc-3 and 
cosines and test to see if the result is of the form 

J'U 1 2 

a. sin (y)Elem(sin (y),cos(y)) 

Substitute 2 = cos (y) 
f 3 ' r 3 

j Sin (y)cos(y)dy becomes j (z - z )dz 

b. cos^'Cy) Elem (cos'* (y), sin(y))dy, substitute Z=sin(y) 

ii) Transform all trigonometric functions into secants and 

tangents and determine if the resulting expression is of 

the form 

a. Elem (tan (y), sec (y)), substitute z * tan(y) 

(sec^ (y) / (1 + sec A (y) - 3 tan(y))dy 


Jsec^ (y)/ (1 + sec a (y) - 3 tan(y))d 
becomes !/(z**- 3z + 2)dz 

nit 1 

b. tan (v) Elem (sec (y) , tan"'(y))dy, 
substitute z * sec (y) 


iii) Substitute z *= tan 1/2 y 


XI. R(x)e where P is a polynomial function of x. See below for a 
description. 


XII. a. In a sum each of the summands is integrated separately. 

b. If the integrand is a small, positive integer power of a :um, 
the integrand is expanded and the result is integrated. 

XIII. (Derivative Divides) See below for description. 

XIV. (Integration by Parts) See below for description. 



Procedure ipt 


r 

Integrating R(x)e 


Algorithm X. 


a “’° mentioned # decision procedure for integrand 
of the foi® K(x)e wnere R la rations, and V is a (non-con*tart) poly, 
nraial. This is an extension of the decision procedure found in Ri't’ -'or 


. T - / \ , , , r - in Kill ;cr 

th^ case I (X; - x. Ths algorithm proceeds by making guesses at the soiut 
Thxs is tne way the Edge heuristic would have proceeded. We ba¬ 
the heuristic slightly to arrive at this procedure 


v e modified 


m, 


Let R(x) *> * where 3^, Q arc polynomials 

Q(x) 5 1 is a polynomial of degree < m 

f 1 «■» M — _ •*» 


is a constant 


** J*» know £ . rco , Rltt that £h * integral (if any) will be a multiple of 
e . Suppose the integral is represented by (^(x) + h L (x))e'’fthen 

+ a^+ •? / (x)-^+ R(x) - cpc "V 


Now, if we guess a t (x) » c' x then 

' P^Q 

/ ft, „ 

a j * nTjC^x - ^ x7 - c, x^P 


P'0 


P' 


and P bj + bj «= R(x) - p'a*- af * Si- m, c, x^+ c, + c ^0 

p' Sa ~ J -TT/-r 


_Pli 


Now, consider the numerator of P ^ + b/ as a polynomial T,(x) 
and a ratronal function remainder U 4 (x). Let the leading term of T( x ) 

<* n Con - lm -. the Process indicated above until come T£ (T r , say) 

0. If at that time the corresponding U,’ (f.e. U.) is also 0 then Me 
expression R(x)e is integrable and the integral is' the sum of the a ‘(xi 
multiplied oy .■«» If U„ U not 0, then the expression is not integribU. 

Note that i. U - 0, then R(x) - p'< a- - is.- =0 
n <• j- i 

Let a » £a; , t.ien we claim that «= 0 implies that 
JR (x ) e^ U) dx * a (> ) e* or 

aP / e^+ a' 

(R - aP'- a')t? ° 

but e'ji 0, henc- R - ap'- a'- 0, but this is assured. 


Example s 
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x-rl 'e*dx is not intejarable 
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Derivative Divides - Algorithm XIII 


Lee the integrand he f(x;g(x) where £(>:) is each of the non- 
constant terms in the integrands considered in turn of the whole 
integrand if it is not a product of terms. 

t f j 

a) If f (x) ® cons tanr*g (x) , then *f(x)g(x) ® f (x) 

' 2*cons cant 

b) If f is a unary operator such a3 sin, cos and its 
argument is k(x) and k / (x) *■ eonstant*g (x) then substitute 
y * k(x) and obtain J 1/constant f(y)dy by looking it up 

in a table. This table contains about ten entries. 

7 

c) If f is a binary operator (i.e., expt or log) then k(x) 
is considered to be the non-constant argument (assuming 
only one argument is non-constant). Proceed as with step 
b above. 

This routine is rather conservative in searching for a subexpression, 
whose derivative divides the integrand. Since it is applied to each 
problem, our main considerations were speed and infallibility in making 
substitutions. Considering its limitations, it is a highly successful 
routine. For example, it helped solve forty-five out of the eighty-six 
SAINT problems with an average time of 0.6 seconds on the 7094. 


Example 


ran r. sec**x dx *= 1/2 tan**x 


s’.** 


e 

AM 


Si a * 


J„ 

Je cos x dx 
jx^dx « l/n-*l x 
Jcos (ax + 1 )dx « 1/e sin(ax *+• b) 



Integration by Parts - Method Xl'v 


At the present time the Edge heuristic discussed in our 
previous memo has net been written. We haue, however, written 
a simple integration by parts routine similar to the one round in 
SAINT. The organization that this routine requires is closer to the 
organization found in SAINT than any other part of the program. 

Consider any partition of the integrand into a product of non¬ 
constant factors g and h, where It * j hdx can be found without call¬ 
ing the integration by parts routine. (This is a modification of the 
SAINT method which required that the integral be found with an allots- 
ment of no resources.) Now consider^Jg»Hdx. We require that this 
integral be found with the number or calls to the integration by parts 
routine fewer than a constant called Maxparts. This constant is 
calculated to be twice the maximum absolute value of an exponent of 
a top level term of the integrand. Thus, xcos(x) may use only two 
calls to integration by parts and x J cos (x) may use four. 

If both integrals can be obtained, then the solution is 
(ghdx * gH - jg'Hdx. 

There are some obvious difficulties with the approach taken 
above. First, there is no attempt to devise a reasonable partition 
of the integrand. Second, there is no routine which determines 
whether progress has been made in obtaining a solution. Both of these 
problems must he solved In a successful implementation of the Edge 
heuristic. 

Note alsc that the value of Maxparts can be set only at the top 
level call to the integration by parts routine. Otherwise, we would 
wind up in an infinite loop. This requires the program to determine 
whether it is at the top level call to this routine. Although it is 
debatable whether we really require such a top-level determination in 
this routine, we predict that complex programs will increasingly make 
use of information regarding the history of the program. For a further 
discussion of this point see Teitelman, 




I; liferent la 1 Equations 


Hourires for solving the following 
equal tons hav.» beer. .written. The ir.tc.rT 
as 2 subroutine in each ccae. 


classes ni differential 
Zlr-r p:u?r;n is called 


a) Separable 

A (x)B (y ) i:< 4- C v .) *> (v} d y * 0 
[a/C dx T yc/p. dy - 0 

x(y*- l)d*. - y(x A - l).vy - o 

becctnc^ ^ >; _ dx - W/ (y*^- l)dv «* C c 

- 1 ^ 

b) Linear 

/+ P(x)y + 0(x) * 0 

becomes ye } •+• lQe'*^dx bj C p 

y / *f y + x ® 0 beccrr.ec vc*t Gce*dx ~ C 

«/ ° 

i) Bernoulli 

y / + P(x)y + Q(x)y 1 * 0 

Substitute v - y* ^ yielding the linear equation (case b) 
v' 4- (l-n)F(x)v 4 (l-n)Q(x) =0 


c) Exa ct 

P(x,y)dx + Q(x,y)dy 
where a Jtfl 


0 

(The matching program SC HAT CHEN is 


o y a x 

used to determine the equivalence.) 
becomes /pdx + Jq - f>dx7)dy « C 

y W 

2xvdx + (x J + cos (y))dy - 0 
becomes x*y + jcos(y)dy = C & 


e) Homogeneous 

P(x,y)dx + Q(x,y)dy * 0 


where P and Q are homogeneous functions of the same degree, 
n, say. Substitute y *» vx and factor out y. n from the 
equation. Substitution yields variables separable (case a). 

(2ye"^ - x)dy -f (2x + y)dx - 0 
y - - 

becomes (2vc - ! )y.dv 4 2 (v’V + 1 )dx - C 



f) Lir/car cocffir.lenis 

y *f F /ex 4- by •*• c_ 

(a'x + b'y + c' 

* s' 


J ' wbere * , c >£,••• #c' are constants 


ar.c a'b - a b'f 0 


substitute x* * x - b'~ 

rrf 


“* D r ' v 

—±-_ ? / 


** *• y - at - a c. 

a'b - 3b' 


a o - an 

and obtain a homogeneous equation (rase e). 


y 


4 y - v + 7 


2x + y - 1 “ 0 bccotnas 4 y*)dy* - (4x - y*)dy* * c 


g) Linear differential equations with constant coefficients 

L. Ernst has recently written a program which solves 
certain cases of inhomogeneous differential equations 
with constant coefficients. The. method used • requires 
Laplace transforms which are obtained by a program 
written by C. Engleman * The latter program uses the 
above-mentioned program by Manove for integrating 
rational functions. We hope to use these programs as 
subroutines. 




Extensions to SCRATCHED 


Two new modes have been added to the m:-. ceding prog can. 
SCHATCHEN. These nave been found us *i : ui in c.ur Lain r> i cua cicn: , 


Consider the expression A{::)E(y>dx 4 Cu)D(v}dy ,{o>.*nd «n 
the section on differential equations when discussing separable 
equations (case a). Vh&t is desired here is that both A and 3 
a product of terms involving x and y respectively. It was not 
previously possible to describe to SCRATCHES a simple matching 
rule which would yield a subset of she terror ir. a product (c.; . f 0 , ; 
that matter„ in a sum). The indicators TIMES** and PLUS** nre'useo- 
for this purpose. For instance* the expression above can be nt-v-ched 
using the following pattern: 

(PLUS (TIMES (TIMES** A FREEY) (TIMES** B FREEX)DX) 

(TIMES (TIMES** C FREEY) (TIMES** D FREEX)DY)) 

where FREEX and FREEY are functions which teat expressions for 
independence of x and y respectively. 


>*> 



rr* 
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REPLACE is to substitute for esch occurrence of u{ in T?;*pv , 
REPLACE will also attempt tc simplify the results obtained ciuri.rt» 
the substitution as the new expression is being formed. R avoids 
certain redundant simplifications b*/ simplifying only at the level 
at which it is working and not at lower levels. (These are assumed 
simplified through the recursion process.) REPLACE also uses two 
special indicators suggestively called Eval* and Ouote*. 

REPLACE(DICX(EVAL* A)) is the simplified result of 

(EVAL(REPLACE DICT A)ALLST). REPLACE(DICT(QUOTE* A)) is simply A. 

By using these indicators we can get some of the dynamic capabilities 
of CONVERT which are lacking in most other string transformation 
languages. 


Example 

(REPLACE 

(QUOTE ((A. - l) (N,l))) 

(QUOTE (PLUS X (TIMES A (EXPT XN)))') ~ 0 
(REPLACE 

(QUOTE ((A, 1) (B. 1) (X.Y)» 

(QUOTE (TIMES X (EVAL* (FACT (PLUS A ?)))))) * (TIKES 2 Y), 

where FaCT is the factorial function. 
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